Vortex Velocity Probability Distributions in Phase Ordering 

Kinetics 

Gene F. Mazenko 
The James Franck Institute and Department of Physics 
The University of Chicago, Chicago, Illinois 60637 
(February 2, 2008) 

Abstract 

The calculation of the point vortex velocity probability distribution func- 
tion (vvpdf) is extended to a larger class of systems beyond the noncon- 
served TDGL model treated earlier. The range is extended to include certain 
anisotropic models and the conserved order parameter case. The vvpdf still 
satisfies scaling with large velocity tails as for the nonconserved isotropic case. 
It is shown that the average vortex speed can be self-consistently expressed 
in terms of correlation functions associated with a Gaussian auxiliary field. 
In the conserved order parameter case the average vortex speed decays as i -1 
compared to the i -1//2 decay for the nonconserved case. 
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I. INTRODUCTION 



In recent work [1] it was shown that the theoretically predicted [2] velocity probability 
distribution for point vortices in the case of a phase ordering system agrees very well with 
direct numerical simulations. In particular the predicted high velocity algebraic tail is found 
to be robust and the predicted exponent confirmed. In the original paper [2] describing the 
theory there were assumptions concerning the Gaussian nature of an underlying auxiliary 
field. Here we clarify this result by showing that the assumption of an underlying Gaussian 
field is consistent and does not imply that the underlying order parameter field is Gaussian. 
We only require that the order parameter field and the auxiliary field share the same zeros 
and symmetry. The theory is also extended here to conserved and anisotropic systems of 
coarsening point defects. 

In ref. [2] it was shown for a nonconserved time-dependent Ginzburg-Landau (TDGL) 
model that for annihilating point defects n = d, where n is the number of components of 
the order parameter and d the spatial dimensionality of the system, that the vortex velocity 
probability distribution function (vvpdf), the probability that a given vortex has the velocity 
V at time t after a quench, is given by 



for long times for the nonconserved TDGL model. Both the form of P(V, t) and the time 
dependence of v(t) have been confirmed in ref. [1] for the case n = d = 2. It is worth pointing 
out that the order parameter growth law, L(t), for the system studied in Ref. [1] has a log 
correction [3,4], L 2 ps t/ln t. This log correction for L(t) is seen in the simulations Ref. [1]. 
There are no log corrections found for v(t). Thus nonlinearities which influence L are not 
seen in v(t). We discuss in more detail here the derivation of the result given by Eq.(l) and 




rXj+i) 
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where the parameter v(t) is clearly related to the average vortex speed and varies as t 1 I 2 
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it's extension to systems with a conserved order parameter and simple spatial anisotropy. 
In the end we find that the vvpdf still satisfies a form similar to Eq.(l) with the same large 
velocity tail, but the average vortex speed falls off as t^ 1 in the conserved case compared to 
the t -1 / 2 behavior found for the nonconserved case. 

The set of problems of interest here are driven by Langevin equations of the form 



dt 

where, we assume, 



KaW (2) 



]im K a ty) = -Oif> a , (3) 

if>— »o 

and the right-hand side is linear in ip. The key idea is that we are interested in that part of the 
equation of motion which corresponds to the motion of vortex cores which are characterized 
by zeros of the order parameter. An important example is the TDGL model which is of the 
form 



dip* = _f 
dt 



V&l>) - cVV« (4) 



where c > 0, f is a constant for a nonconserved order parameter (NCOP) and f = —DV 2 
for a conserved order parameter (COP). Comparing Eq.(4) and Eq.(3) we have 

d NCO p(i) = -r c v 2 (5) 

and 

OcopiX) = DV\[-r + cV 2 ] (6) 

for the COP case. Here r = V"(t/j)\^, =0 and r < if the system is unstable. Through a 
proper choice of length and time scales we can choose 

Oncop(1) = -V 2 (7) 



<W(1) = V?[1 + V?] . (8) 

An anisotropic model can be assumed to be of the form 

6 ANI = -£ c M vj, + £ ^,m 2 vj, Vj 2 (9) 

where c M1 and & Ml)M2 are constants. This reduces to the COP case if c M1 = — 1 and & MliM2 = 1- 
There is the underlying assumption that the nonlinear potential contribution in the 
equation of motion must be such that system orders via annihilating point defects. 

II. DEFECT DENSITIES AND CONTINUITY EQUATION 

We assume that the instantaneous positions of these defects are determined by the zeros 
of the order parameter field. Furthermore, it was pointed out in ref. [2], that the vortex 
charge density for this system can be written as p = 5(ip)V where V is the Jacobian (deter- 
minant) for the change of variables from the set of vortex positions ri(t) (where vanishes) 

— # 

to the field ip: 

^ = _ t e ^l,V2,---,Vn e Vl,V2,---,Vn^ Vl^Vl^ V2^V2----'V fln^Vn (10) 

Til . 

where e Ml)M2) ... )/ln is the n-dimensional fully anti-symmetric tensor and summation over re- 
peated indices here and below is implied. Furthermore, since topological charge is conserved, 
it was shown in ref. [2] that p satisfies a continuity equation: 

| = -V.(pv) (11) 

where the vortex velocity is given by 

VV/3 = ~ ( n _ 1 ^| e ^2,.,Mn^l^2,...,^(O^l)V )U2 -0^----V Mrt -0^ • ( 12 ) 

It is assumed that the velocity field is used inside expressions multiplied by the vortex 
locating 5-function so we can use Eq.(3) in Eq.(12). These results are rather general. 
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Notice that p and v have certain important invariance properties. If we can write 

^„(1) = (a + (3m\l) + ...)m u (l) (13) 

for small m v and where a and (3 are constants, it is easy to see that p(ip) = p(rh). If we 
further assume, as m and ip go to zero, 

6^(l) = a6m„(l) , (14) 

then 

— * 

= v^m) . (15) 

If corresponds to an operator with two gradients, as in the NCOP TDGL model, then, as- 
suming Eq.(13) to be valid, Eq.(14) follows. In the case where O has higher-order derivatives 
and we assume Eq.(14) holds, then Eq.(13) must be modified. 

Thus the correlation function we compute in the next section, (2(12), is for that set of 
fields m, related to to ip by Eq.(13), for small values of m v and ip v , which is described by 
a Gaussian distribution. Thus we assume there is a field m v which is Gaussian while the 
statistics of ip u are largely undetermined. 

III. CORRELATIONS FOR THE DEFECT SECTOR 

We show in Sect. IV that in determining the average vortex speed we need certain 
correlation functions for the auxiliary field, m, introduced in the last section. We show here 
that we can use the defect continuity equation, Eq.(ll), to show that there is a self-consistent 
solution for an m field that is Gaussian. Furthermore we determine the correlation functions 
needed to evaluate the average vortex speed explicitly. The method we develop here is a 
generalization of the approach due to Mazenko and Wickham [5] . 

The idea is to look at the equation generated by multiplying the continuity equation by 
a source function and then averaging over m. We have 
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([dp(l)dt l + vW-(p(lMl))}s(H)) = (16) 

where S(H) = exp[J rflH(l) ■ m(I)]. The question is whether this equation can be satisfied 
for an underlying Gaussian probability distribution for arbitrary external field H(l)? 
To answer this question we evaluate first the quantity 

(p(l)S(H)) = (S(H)8(1)V(1)) (17) 

where we introduce the simplifying notation 6(1) = 5(m(l)) and now T> is a function of 
the field m. When we talk about correlation in the defect sector we mean averages like in 
Eq.(17) where there is a vortex locating 5-function inside the average. 

By taking functional derivatives we are able to generate the correlations between fields 
m at arbitrary space-time points with the field at the space-time point 1. If we define 

Z H (1) = (S(H)S(l)) (18) 

then 

(^(1)^(2)^(3) . . .) = Z- H ^) m S n(2) fg * (8) . . . Z„(l) . (19) 
In our development a key property of the underlying Gaussian distribution function is 

for arbitrary A. For A = m U2 {2) we obtain 

(m Vl {l)m V2 {2)) = G VlV2 (\2) . (20) 

If we assume that the system is isotropic in the vector space, then G UlU2 (12) = 5 VlV2 G(12) 
and 

{mnim= j iiam^A) . (2D 



We then need to work out 

(S(H)p(l)) = (S(H)S(1) ) • (22) 

Using Eq.(21) we have 
(S(H)p(l)) = 

G(l!)( ^(T) (5(iJ)5(1)v ' t2m ^---- v ^ m ^ )) • (23) 

The derivative of S(H) leads to the introduction of the quantity 

A U1 (1)= j dlG{ll)H vl {\) . (24) 

We assume, and check self-consistently, that (V Atl G(ll)) |i = i = 0. The derivatives of the 
product V ^ 2 m V2 ....V ^ n m Vn with respect to m Ul (I) lead to contributions which all vanish 
because it picks out terms S vlVj which multiplies e Ul ,v2,.,vj,.,v„ an d e vi,v2,.,v u .,v n — 0- We have 
then 

(S(H)p(l)) = V ltl A in (l){S(H)6(l)V lta m va ....V lin m Vn ) . (25) 

Clearly we can go through this process n — 1 more times to obtain (S(H)p(l)) = T>a(1)Z h (1) 
where 

V A {1) = ^eMi,M 2 ,...^^ 1 ,, 2 ,...,, n V w ^ 1 (l)V M2 A ;/2 (l)....V M „^ n (l) . (26) 
Next we look at the current contributions to Eq.(16) in the form 
(S(tf)V>(l)v(l))) = -V«--— (S{H)5{l)fn v ,V» 2 m V2 ....V lin m Vn ) . (27) 

We assume that 

5(1)771^(1) = -5(1)0(1^(1) (28) 

where 0(1) is a derivative operator defined by Eq.(3). In ref. [2] one has the choice 0(1) = 
— Vf. After inserting Eq.(28) back into Eq.(27) and using Eq.(21) we obtain 
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(S(H)V- (p(l)v(l))) = -VS^^,^ ,,^, 2 % 



x 



(-6(1)^(1)) V, 2 A V2 {l)....V^A Vn {l)Z H {l) 



+ (-6(l)G(ll)) i= _ (5(iJ)^ 1 (l)V M2 m i , 2 ....V M „m i ,J 



(29) 



where 5 V1 (1) = dm d ^ 5{m(l)) and the derivatives of the product of m's vanish as in the pre- 
vious case. Clearly the reduction of the term containing 5 Vl (l) follows just as for (S(H)p(l)) 
with the result: 

(S(H)V ■ (p(l)v(l))) = V^^^e^,...,^^,...^ 



x 



(0(1)^(1)) V lta A V2 (l)....V lln A Vn (l)Z H (l) 



+ (6(l)G(ll)) i= _V w ^(l)....V /lB ^(l)(5(^)^(l)> 



(30) 



We must next evaluate Zh(1) and the related quantity (S(H)S V1 (1)). Determination of 
ZuiX) involves evaluation of the Gaussian average 

d d k 



Z H (1) = (<5(m(l))S(H)) = / ^J^^Mi) e m)Mi)) 



with the result 



e 2 s oW 
z h(1) = - - ,^„ n/9 exp 



(2nS (l)) n/2 
where S (l) = G(ll). Next we need 



(2?r) f 



-^(1)^(2)^(12) 



d d k 



(S(H)5 U1 (1)) = J i^^ 1 )^)^ 1 )) 



(31) 



(32) 



d d k 



(2k) 



~k 2 G(ll) + *k ■ A(l) + l -H V2 {l)H V2 {2)G{\2) 
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exp 



-H Vl {l)H Vl {2)G{\2) 



d 



dA Vl {l) 



exp 



--H Ul {l)H Vl {2)G{\2) 



Zh(1) 



AM 



Zh(1) 



Putting Eqs.(32) and (33) back into Eq.(30), and allowing the gradient to act 
(S(H)V ■ (p(l)v(l))) = -V B (1)Z H (1) -V B (1)V%Z H (1) 



where 



Vb{1) 



(n-1) 



■£(ii,/i 2 ,...,fi„ e i'i,i'2,...,i' n V w 5 i / 1 (1) V M2 A U2 {1) . . . .'SI ^ n A Vn {l) 



and 



i 

(n- 1)! 



s^(i) = (-d(i) + n(i)) 



and 



= 5^1) (Ot^ 12 )). 



Putting the results together in Eq.(16) we obtain 



J- (^(1)^(1)) = 2> fl (l)Z H (l) + ^(l)V«Z ff (l) . 



We can write this in the form: 



After taking the derivatives Z#(l), we can then write Eq.(40) in the form 



W 2 (l) = W 4 (l) 



where 



and 



Look first at W2(l) which can be written in the form: 



1 

(n-1)! 



ffii.ft,..,/!,^!,^,..,^ ^ m9vi iX) V^-A^ (!)•••• V Mn Aj, n (l) 



(n-1) 



where 



^(l)=^(l)- ^|^y^(l) -^(1) 

In looking at W 4 we need to focus on the quantity 

2£(1)4,(1)V W 4,(1) 

1 



~ (n-1)! 
Note that 



S J/1 (l)V M2 ^ 2 (l)....V Mn ^ ?i (l)^(l)V m ^(l) 



e w ^,...,^V m ^(l)V^^ 2 (l)....V Mn ^ n (l) = ,,'P.i(U . 

Putting this back into Eq.(46), we find 

^ i (l)^(l)V w ^(l) = ^ I ^e, 1 ,, 2 ,..., i/n e^ 2 ,...,,„ J B i , 1 (1)^,(1)^(1) 
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= B Vl {l)A Vl {l)V A {l) 



(48) 



and 

= V ^§^sM ■ (49) 

Using Eqs.(44) and (49) in Eq.(41) we find a solution for general H if 

^(l)=i,(l)-i|^|^(l)-^(l) = . (50) 

This will hold for all source fields H if 

^ G ( 12 )-\^ G ( 12 ) = (-0(l) + m)G(l2) . (51) 

Thus the average of the continuity equation, Eq.(16), is satisfied by a Gaussian probability 
distribution if the associated variance, G(12), satisfies Eq.(51). 



IV. SOLUTION FOR G(12) 

We can solve Eq.(51) for G(12) in some generality. The first step is to write 

G(12) = v /5o(l)5 (2)/(12) (52) 

where So(l) = G(ll). Inserting this form into Eq.(51) gives 

A /( i2) = (-(5(l) + n(l))/(12) . (53) 

We assume that the system is translationally invariant and on Fourier transformation the 
operator 0(1) is diagonalized and time independent. We have then 

-^f(q,t 1 ,t 2 ) = (-0(q) + n(l))f(q,t 1 ,t 2 ) . (54) 
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We see that Q(l) , defined by Eq.(38), is determined by the constraint 

d d q 



m = J J^ d O{q)f{qMM) ■ (55) 
We also have the equation 

J?-f(q,t 1 ,t 2 ) = (-0(q) + n(2))f(q,t 1 ,t 2 ) . (56) 
ot 2 

Adding Eqs.(54) and (56) and setting t\ — t 2 — t we obtain for the equal-time correlation 
function f(q,t) = f(q,t,t): 

p(q,t)=2(-0(q) + n(2))f(q,t) . (57) 



For equal times we have from Eq.(52) that /(ll) = 1 or 

1 = I ■ (58) 

The partial solution for Eq.(57) is given by 

f(q, t) = exp (2 £ dr (Q(r) - 0(g))) f(q, t ) 

= R 2 (t,t )e- 2O ^ t ^f(q,t ) (59) 

where 

R(t u t 2 ) = exp (jf** dr Q(r)J . (60) 
We then need to solve for Q(t). Inserting Eq.(59) into Eq.(58) gives 

1 = R 2 (t,t )I(t,t ) (61) 

where 

/(Mo) = / ^ e ' 2 ° i9){t ~ t0) f^ ■ ^ 



d d q _ 2 
(2^) c 
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We then have from Eq.(61) 

R\t,t ) =I~\t,t ) . (63) 
The constraint condition, Eq.(55), is given by 

m=R 2 (t,t ) J ^ d O( q )e-^-^f( q ,t ) 



ii(t,t ) 



(64) 



2 /(Mo) 

Thus the determination of Q(l) is reduced to evaluation of the integral I(t,t ). The equal 
time correlation function is given then by 

f(q,t) = r\t,t )e- 2O ^-^f(q,t ) . (65) 

Going back to the unequal time correlation function we can integrate Eq.(54). 

f(q, h, t 2 ) = exp dr (Q(r) - 0(g))) f(q, t 2 ) 

= R(t u t 2 )e-°^'^f(q,t 2 ) 

= R(h, t Q )R(t 2 , t )e'°^ +t ^f(q, t ) (66) 

where we have used R(ti,t 2 ) = R(ti,t ) / R(t 2 ,t ) . We obtain a complete solution once one 
does the integral for I(t,t ). Notice that these results are independent of the specific form 
for S (t). 

V. EVALUATION OF VORTEX VELOCITY PROBABILITY DISTRIBUTION 

The vortex velocity probability distribution function defined by 

noP(V,*) = (n(l)<J(V-v(l))) (67) 
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where v as a function of the order parameter is given by Eq.(12) with if replaced by to, 
n(l) = 5(m(l))|P(m(l))| is the unsigned defect density, and n (l) = (n(l)). We notice 
in evaluating P(V,t) that it is of the defect sector form, thus there is a defect locating 
(^-function in the average via the factor of n(l). Our results from section 2 suggest that in 
this sector we can treat the field m as Gaussian with variance given by (2(12) calculated in 
section 3. 

In carrying out the average we need the auxiliary quantity 



b) = (5(m) J] - V,m u )5(h - K) (68) 

fl,U 

where K(l) = 0(l)m(l). Then 

n P(V,t) = [ d n bl[dZ;\V(Z)\5(V-v(b,Z))W(Z,b) (69) 



where v(b, 0) = J(b, 0/^(0 with 



"^(0 — i^l^j,..,)!,^!.^,../,^^""^ (^) 



ni 
and 

Let us turn to the Gaussian average determining W(£,b). Following the analysis used 
to evaluate Z H in section 2 we introduce the integral representation for the 5-function to 
obtain: 



where 

r(k,q,s) = (e lk - m(1) e iq - K(1) e ls M V M"^(i)) _ ^ 

If we introduce 
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H a {\) =i[k a + q a 6(l) + s^V«] 5(11) 

then we can write 

r(k,q,s) = ( e ^(i)™ Q (i)) = exp 



l -H a {l)H a {2)G{Y2) 



where G{Y1) was determined in section 3. We assume that the cross terms involving 
number of gradients vanishes in the argument of the exponential. We have then 



r(k, q, s) = exp 

where 



-i {k'S a (i) + 2k ■ qSyi) + , 2 S 02 (1) + s«s«Sf>(l 



50(1) J (2tt)« wj v *' ' w 2/(1) ' 
where we have used Eq.(54), 

5 02 (1) _ f d n q ... . _ 1/(1) 



5 (1) / (2tt)" y 4/(1) 

having used Eq.(65), and 



The next step in extracting is to integrate over k: 



where 

r< o c 

<J — >->02 — 77" • 

Using Eq.(80) back in Eq.(72) we obtain 
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This factorizes into a product of three natural parts 

^' b)= (^p W( ^ (b) (83) 



where 



W(0 = (il/ ^f) e-«e-^4 2) (84) 



and 



Using the basic integral 



^ (b) = / Wr e ~ lh '" e ~ h2 ' s ■ (85) 



^e-^e-f- 2 = _L e "£ (86) 



we can evaluate both factors: 



/ \ n/2 _ " 

^(n^) e ^ ■ < 87 > 

H/(b) = (S^ 6 • (88 » 

Turning to n P(V,t) given by Eq.(69), we see that we have the integral over b of the 
form 



J b = J d n b5(V-v(b,0)W(h) 



I dnh I ?fi e ~ iV ' ZeMb ' ' zw ( b ) • ( 89 ) 



(2w) n 

We can then write z • v(b, £) = a ■ b where 
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X>(n — iy Zc * ea,ll2 '---^ nCui '' /2 '--- , ' /n ^L l2 ""^L ln 



£V2 tv n 



(90) 



then 



J, 



(2vr)« 



e~ iV ' z e ia ' b 



1 



(2nS) n l 2 



e as 



r d n z 
~ 1 (27r) r 



e -,V.z e -fa2 



where a 2 = z a M a pZp and the matrix M is given by 
1 



v 2 [(n- ly.y 



, , t^2 C V ^(z , ,z , t"2 t^rj 

^2v",^ C W,...,l'nS / i 2 '-'-V/^ 2 ,.,<V2v",l'! 1 < >|l' ""V!, 



Doing the remaining Gaussian z integration in Eq.(91) we obtain 



and 



(91) 



(92) 



(93) 



(94) 



We must look at the matrix M and its inverse. First multiply M a/3 by ^ to obtain 
1 



However 



V 2 [{n- 1)!] 



rea,/i2v--^n e ^l,^2v--^n^^---^™ e /3,^,...,^ e ^l,^2v--,<^----^ ' ' 



(96) 



then 



i a M a p — jy 2 ^ n _ ]_^|j2^^) e ^^ 2 '---'^ e ^ 1 ^ 2 '---^™ e ^'A t 2v--^n e ^1^2v--:<^4''''^n 



(97) 



(98) 
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where we have used 



^v,V2,...,v„^vi,i^i,...,v n 3i/,i>i{ n 1)' • (99) 

Multiply Eq.(98) by ^ to obtain 

= V ( n _ ^ ^...x^tO^ i<, = 5 7/3 • (ioo) 

Thus we have the beautiful result: 

{m-% = E^ • ( 101 ) 



We need det M = 1/det M~ l . We have 



detM 1 — —^ ai ,a2,...,a n tl3 1 ,P2,...,PnCaA7A V a 2 Q 2 2 ' ' ' ZaZQl 



= ^V(0 2 n\=V(0 2 ■ (102) 



Using the clean result det(M) = 1/(V) 2 , and Eq.(lOl), we have 

Vr d% \ V\Q 



where 



^) = Ei(C) 2 4Em/ . (104) 



Next make the change of variables f£ = V 'S^Q to obtain, £>(£) = V^ 2) ) £>(£)> 
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n P(V,t)=n 



S (2) 



n 



U(0 



where 



Mi) = E© 2 + E v a tiW 

W a,f3,v 



and 



V c 



c(2) 
Oct 



V a . 



Next we make the transformation from to x« via 



(1 ! 



la = 



(I' 



such that 



Mi) = EM 2 = E^^&x* + E ^^^x^ 



This requires that iV satisfy 



(I 1 



,/*2 "Mil. 



M2 



(1 



This has a solution given by 



We then need the Jacobian of the transformation 
terms of x- Look first at 



c?X^ and X>(£) evaluated 



^(0 ~~ ■ Eai«2...a„E|/ii/2...f n 
?7>. 



X 
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If we multiply this out in powers of V we see that if we have more than one factor of V a% 
then the contribution vanishes due to antisymmetry, therefore 



X>(£) = V(x) + - } TKa iaa ..^w..^N 1 V™rf 1 V*x£ . . . xZ 



which can be put in the form: 



m = 



gfr) 

\Jl + V 2 



Next we need the Jacobian 



J = Y[ det 



K 



(113) 



(114) 



(115) 



det (S^ + N^V^Y = (J ) 



(116) 



where 



Jo = det + 



(117) 



1 

n\ 



| e aiQ2---Qn e /3l/32.../3n 



S aiA + Nrf^V*] \S aatfh + N x V°*Vfr] . . . \5 a:iA + N^V*} . (118) 



Again, expanding this out in powers of V, only the first two terms contribute due to sym- 
metry and we have 

Jo = ^/ ai a 2 ... an + nN 1 e aia2 ... an e Pia2 ... a S ai V 131 



1 + 7Vi\/ 2 



Going back to Eq.(105) we have 

S {2) 

n P[V,t] = l[ l " 



l + \/ 2 



,\2n^ (l + V*)(»+2)/2 



J* 



(119) 



(120) 
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where we have the final integral 



\Mx) 



121) 



We can evaluate Jf directly. The first step is to write: 



Jh 



/n 



e e>, , y (1) y (2) y {n) y {x ) y {2 ) y^e"^,^) 2 

t^l/i 2 ...Atn t M>2---^nA./il ^ H2 ' ' ' A,l„ A^ A^ ■ ■ ■ Xfl' n C 



122) 



This factorizes into a product of integrals for fixed v 

dxW m (i) 4V fv^^ /•i-r rf xL 2) 



y y ^ ^ x ^ y y v^f w K 



Each integral in the product is equal to 1 except for those giving a 5-function with unit 
coefficient: 



Jp = e 



— <ri I 

^7UM2— /in 



(123) 



We have then 



n P[V,t] = n!]J ! " 



V V27TVW/ (l+^ 2 )(" +2 )/ 2 

Since P(V, t) is normalized to one, we find on integration over V, the result 



(124) 



n 



ni 



2 n / 2 Y{ 



o(2) 
Ofj, 



(125) 



which agrees with previous results in the isotropic limit. Eliminating n Q in Eq.(124) we 
obtain the final result for the vvpdf: 

1 



7T 



^Vj (1 + \/ 2 )(™+ 2 )/ 2 



(126) 
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where V a = V a /v a , 



using Eqs.(77), (78) and (81) we have 

1 



S = S - 



while Eq.(79) gives 

with I defined by Eq.(51): 



(127) 



(128) 



(129) 



(130) 



The needed input to determine the average vortex speed, v^, is the function O(q), the 
Fourier transform of the operator defined by Eq.(3), and the initial condition /(q, to). This 
result for P[V, t] is the anisotropic generalization of Eq.(l). For the set of models included 
here the velocity tail exponent is (n + 2) independent of direction. 



VI. ANISOTROPIC CASE 

As a particularly simple example, suppose that we have the choice in the governing 
Langevin equation 

0(l)^(l) = -EcaV^(l) , (131) 

a 

or, in terms of Fourier transforms, 

0(q)=J2c a ql ■ (132) 

a 

The associated vortex- velocity probability distribution is given by Eq.(126) and we need to 
work out the average vortex speed given by Eq.(127). Assuming the initial condition 
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/(g,0)= (n(2^ Q ) 1/2 )e-M , 

we have from Eq.(130) 
where 

h a {t) = h a + Ac a (t-t ) , 

from Eq.(128) 
and from Eq.(129) 

f d n q 2x/ _, 1 



S (l) J (2tt)"^ w A ' h,it) 

Putting these results back into Eq.(127) we find the scaling velocity for a 
system is given by 



vl = 2Kit)J2 



\ 2 



f VM*). 



( (*-t ). 



In the large time limit we have 



_ 2 d 
vf, = — - 
2 t 



and the final form is a simple generalization of the isotropic result. 
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VII. CONSERVED ORDER PARAMETER CASE 



Let us look at the COP case where 0(q) = —q 2 + q A . We choose the rather general initial 
condition 

/ h \ d/2 

/M)=(^-J e-^ 2 , (141) 

which satisfies the normalization condition given by Eq.(58). We then need to evaluate the 
integral I and the numerator in Eq.(129): 

r d d q q 2 2t(g2 _ g4) (h_\ d/2 ^ 
J (2n) d d \2n) 

= -h d / 2 ^r(h- d / 2 I) . (142) 

(Jj CJ I V 

We see that all of the ingredients contributing to v 2 can be expressed in terms of I and it's 
derivatives. We see that / can be written in the form: 

/ = I Q h d ' 2 r dqq^e 2 ^-^ 1 --^ (143) 
Jo 

where Jo is a constant which depends on d and cancels when we take ratios. Changing 
integration variables to x = q 2 we find 

/ = i'tfl* f°° dxxW-W*-^*-** (144) 
J o 

where I' Q = Iq/2. The leading large time dependence can be extracted from the integral 
by completing the square in the argument of the exponential or using the stationary phase 
method. We find, to leading order in large t: 

' , \ d/2-l . 



/ = / °UJ Va** (145) 

where b = 1 — ^, 2 = |6 2 . From this result for / we see that I = ujI where 

« = <i-i)s + 2tf-s = 5-s + 0(rl > ' (146) 
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Going further we have I = uj 2 I + Col which leads easily to the useful result: 



I 
1 



(147) 



Turning to the evaluation of J given by Eq.(142), using Eq.(145), we find 



J = — 

2dt 



(d/2 - 1) 



(148) 



Working to leading order in time we find J = ^. Putting all of this together in Eq.(127) 

9 Co d . 

11 = m = 7- (149) 

since, from Eq.(146), Co = + 0(t~ 3 ). The final result for v 2 is independent of the initial 
conditions. We see that the COP average vortex speed is qualitatively slower than the 
NCOP case: 



v cop _ 1 



Vncop t 

The computation of using Eq.(127) has been checked numerically in the simplest n 
d = 2 case where / can be evaluated explicitly in terms of an erfc function. 



(150) 



VIII. CONCLUSIONS 

We have presented here the detailed calculation of the vvpdf including the time de- 
pendent vortex scaling velocity for a class of models beyond the original nonconserved 
TDGL models. The class of models studied includes the conserved TDGL model and certain 
anisotropic models. In the conserved case it is found that the average vortex speed falls off 
as t^ 1 compared to the NCOP case where v ~ t~ 1 / 2 . It is our intension to numerically test 
the predictions for the COP case. 

We see that there is self-consistent confirmation that in dealing with vortex velocities 
one can organize things in terms of averages over an auxiliary Gaussian field. We require 
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self- consistently that this field and the order parameter field share the same zeros. A similar 
development can be worked out for string defects [5,6]. 
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